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We present measurements of the branching fractions and charge asymmetries (where appropriate) 
of two-body B decays to rf'^K* , rf'^ p, n^n , tun , and (fin . The data were recorded with the BABAR 
detector at PEP-II and correspond to 89 x 10 BB pairs produced in e + e~ annihilation through 
the T(4S) resonance. We find significant signals for two decay modes and measure the branching 
fractions B(B+ -> r)K* + ) = (25.6±4.0±2.4 ) x 1(T 6 and B(B° -> nK* ) = (18.6±2.3± 1.2) x 1(T 6 , 
where the first error is statistical and the second systematic. We also find evidence with significance 
3.5<t for a third decay mode and measure B(B + — > r/p + ) = (9.2±3.4±1.0) x 10~ 6 . For other channels, 



G 



we set 90% C.L. upper limits of B(B° -> r/p°) < 1.5 x 10 -6 , B(B+ -> ■q'K* + ) < 14 x 10~ 6 , B{B° -> 

n'K*°) < 7.6 x icr 8 , b(b+ -> V ' P + ) < 22 x irr 6 , b{b° -> v 'p°) < 4.3 x nr 6 , b(b° -> 7771-°) < 

2.5 x 10" 6 , B(B° -> r/'n ) < 3.7 x 1(T 6 , B(B° -> wtt ) < 1.2 x 10" 6 , and B(B° -> 07T°) < 1.0 x 10~ 6 . 
For self-flavor-tagging modes with significant signals, the time-integrated charge asymmetries are 
A ch (riK* + ) = +0.13 ± 0.14 ± 0.02 and A ch (r]K*") = +0.02 ± 0.11 ± 0.02. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

We report the results of searches for charged or neutral 
-B-meson decays to the charmless final states [1] ijK*, 
r]'K* , T]p, r/'p, 777T , ij'n , ujtt , and ^>7T°. For decays that 
are self-tagging with respect to the b or b flavor, we also 
measure the direct CP-violating time-integrated charge 
asymmetry, 



Ach = 



r- -r+ 
r- + r+ 



(i) 



The superscript on T corresponds to the sign of the B ± 
meson or the sign of the charged kaon for B° decays. 
Throughout this paper, we use 77 W to indicate either r\ 
or rf . 

Interest in B decays to r\ or 77' final states intensified in 
1997 with the CLEO observation of the decay B -> r/'K 
[2]. It had been pointed out by Lipkin six years earlier 
[3] that interference between two penguin diagrams (see 



W+ 



B+ u,c,t 




B+ u,c,t 




Figs, la and lb) and the known 77/77' mixing angle con- 
spire to greatly enhance B — > rj' K and suppress B — > r\K . 
Because the vector K* has the opposite parity from the 
kaon, the situation is reversed for the B — ► r/'K* and 
B — * 77-RT* decays. The general features of this picture 
have already been verified by previous measurements and 
limits. However, the details and possible contribution of 
the flavor-singlet diagram (Fig. Id) can only be tested 
with the measurement of the branching fractions of all 
four (r/,ri')(K, K*) decays; the branching fraction of the 
B — > rj'K* decay is expected to be particularly sensi- 
tive to a flavor-singlet component [4, 5]. The tree dia- 
gram (Fig. lc) is suppressed by the parameter A of the 
Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix. 

By contrast, for the B — > rj^p + decays, the penguin 
diagrams (Figs. 2c and 2d) are CKM- suppressed. Since 
the internal tree diagram (Fig. 2b) is color-suppressed, 
the decay is dominated by the (external) tree diagram of 
Fig. 2a. 

The B° decays are different because there are no ex- 
ternal tree diagrams analogous to Fig. 2a. In Figs. 3a 
and 3b we show the penguin diagrams and in Figs. 3c 
and 3d the color-suppressed tree diagrams for the B° — > 
T1 (')p0 : b° -> 77071-°, and B n -> wvr decays. The 
color-suppressed diagrams cancel for the rj and 77' de- 
cays and are expected to be largely suppressed for the 



pseudoscalar-vector (PV) B° 



decay. The singlet 



77,77 



penguin diagram (Fig. 3e) may be significant only for 



(a) 



(b) 




(d) 



FIG. 1: Feynman diagrams for the decays Bp —> (r],T]')K* + . 
The corresponding neutral decays are similar except that the 
spectator quark becomes a d, the gluon in (b) makes dd, and 
the tree diagram in (c) has an internal W. 



*Also with Universita della Basilicata, Potenza, Italy 
tAlso with IFIC, Instituto de Fi'sica Corpuscular, 
Universidad de Valencia, Valencia, Spain 
t Deceased 



CSIC- 



(c) 



(d) 




FIG. 2: Feynman diagrams for the decays B + — » r/p + and 
B+^r,'p+. 
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the decays with an rj in the final state, and the elec- 
troweak penguin (Fig. 3f) is the only contribution for 
the B° — > 4>ir° decay (and negligible for the other decay 
modes) . Branching fractions for all these decays are gen- 
erally expected to be in the range (0.1-10) xlO -6 [6-9], 
with the B + — ► r/^ p + decays at the high end of this 
range and the B° decays at the low end (and B° — > ipn 
perhaps somewhat below this range). 

The charge asymmetry A c h for most of these decays is 
expected to be < 10% [6, 10]. However, for B -> rj ' K* 
the penguin and tree amplitudes are expected to be 
of similar magnitude, which allows charge asymmetries 
which could be in the (20 - 40)% range [5, 8, 9, 11]. 
Information on charge asymmetries and branching frac- 
tions from this full collection of B decays can serve to 
constrain the relationship between the various underly- 
ing amplitudes. 

The results described in this paper complete the mea- 
surement of all four (rj,ri')(K, K*) final states, as well 
as those with (77, r)')(TT,p), with a BABAR dataset of 89 
million BB decays. Current knowledge of the decays 
discussed here comes from published measurements from 
CLEO [12-14] and BABAR [15]. Results for the final 
states (rj, r]')(K, w) on this dataset have been presented 
elsewhere [16, 17]. These data represent an order of mag- 
nitude increase in the B meson sample size over the only 
previous complete study. 

All results are based on extended maximum likelihood 
(ML) fits as described in Section V. In each analysis, 
loose criteria are used to select events likely to contain 
the desired signal B decay. A fit to kinematic and topo- 



W+ W+ 




b u b u 




d d d d 




00 (f) 

FIG. 3: Feynman diagrams for the B° decays. 



logical discriminating variables is used to differentiate 
between signal and background events and to determine 
signal event yields and time-integrated rate asymmetries. 
In all of the decays analyzed, the background is dom- 
inated by random particle combinations in continuum 
(e + e~ — > qq 1 q = u,d, s,c ) events. Some decay modes 
also suffer backgrounds from other charmless B decays 
with topologies similar to that of the signal. In such 
cases, these backgrounds are accounted for explicitly in 
the fit as discussed in Sec. IV C. Signal event yields are 
converted into branching fractions via selection efficien- 
cies determined from Monte Carlo simulations of the sig- 
nal as well as auxiliary studies of the data. The complete 
analysis is carried out without regard to whether there 
are observed signals. This "blind" procedure is used to 
avoid bias in the results. 



II. DETECTOR AND DATA 

The results presented in this paper are based on data 
collected with the BABAR detector [18] at the PEP-II 
asymmetric-energy e + e~ collider [19] located at the Stan- 
ford Linear Accelerator Center. The results in this pa- 
per correspond to an accumulated integrated luminosity 
of approximately 82 fb _1 , corresponding to 89 million 
BB pairs, recorded at the T(AS) resonance ("on-peak", 
center-of-mass energy y/s = 10.58 GeV). An additional 
9.6 fb _1 were recorded about 40 MeV below this energy 
( "off-peak" ) for the study of continuum backgrounds in 
which a light or charm quark pair is produced. 

The asymmetric beam configuration in the laboratory 
frame provides a boost of f3-f = 0.56 to the T(4S). 
This results in a charged-particle laboratory momen- 
tum spectrum from B decays with an endpoint near 
4 GeV. Charged particles are detected and their mo- 
menta measured by the combination of a silicon vertex 
tracker (SVT), consisting of five layers of double-sided 
detectors, and a 40-layer central drift chamber, both op- 
erating in the 1.5-T magnetic field of a solenoid. The 
transverse momentum resolution for the combined track- 
ing system is <t Pt /pt — 0.0013pT © 0.0045, where the 
sum is in quadrature and px is measured in GeV. For 
charged particles within the detector acceptance result- 
ing from the B decays studied in this paper, the aver- 
age detection efficiency is in excess of 96% per particle. 
Photons are detected and their energies measured by a 
CsI(Tl) electromagnetic calorimeter (EMC). The photon 
energy resolution is a E /E = {2.3/ E( GeV) 1 / 4 1.9} %, 
and the angular resolution from the interaction point is 
ag = 3.9°/ \JE{ GeV). The photon energy scale is deter- 
mined using symmetric 7r° — > 77 decays. The measured 
7T° mass resolution for 7r°'s with laboratory momentum 
in excess of 1 GeV is approximately 8 MeV. 

Charged-particle identification (PID) is provided by 
the average energy loss (dE/dx) in the tracking devices 
and by an internally reflecting ring-imaging Cherenkov 
detector (DIRC) covering the central region. The dE/dx 
resolution from the drift chamber is typically about 7.5% 



8 



for pions. The Cherenkov angle resolution of the DIRC 
is measured to be 2.4 mrad, which provides a nearly 3tr 
separation between charged kaons and pions at a mo- 
mentum of 3 GeV. Additional information that we use 
to identify and reject electrons and muons is provided by 
the EMC and the detectors of the solenoid flux return 
(IFR). 



III. CANDIDATE RECONSTRUCTION AND B 
MESON SELECTION 



u>, and selection 



We reconstruct B mesons in the final states r]^K* + , 
rj^K* , r/')p+, r/V, t/^tt , W7T°, and ^tt°. Monte 
Carlo (MC) simulations [20] of the signal decay modes 
and of continuum and BB backgrounds, and data control 
samples of similar modes, are used to establish the event 
selection criteria. The selection is designed to achieve 
high efficiency and retain sufficient sidebands in the dis- 
criminating variables to characterize the background for 
subsequent fitting. As the invariant mass distributions 
from the primary resonances (rj^ , K*, p, uj, and 4>) in 
the decay are included in the maximum likelihood fit, the 
selection criteria are generally loose. Additional states — 
7T° or 77 in f]' decays, and K® — are selected with the re- 
quirement that the invariant mass lie within 2-3tr of the 
known mass. 



A. Charged track selection 

We require all charged-particle tracks (except for those 
from the — > 7r+7r~ decay) used in reconstructing the 
B candidate to include at least twelve point measure- 
ments in the drift chamber, lie in the polar angle range 
0.41 < 8iab < 2.54 rad, and originate from within 1.5 
cm in the x — y plane and 10 cm in the z direction from 
the nominal beam spot. We require the tracks to have a 
transverse momentum pr of at least 100 MeV. 

We also place requirements on particle identification 
criteria. We veto leptons from our samples by demand- 
ing that tracks have DIRC, EMC and IFR signatures 
that are inconsistent with either electrons or muons. The 
remaining tracks are assigned as either charged pion or 
kaon candidates. This assignment is based on a like- 
lihood selection developed from dE /Ax and Cherenkov 
angle information from the tracking detectors and DIRC, 
respectively. For the typical laboratory momentum spec- 
trum of the signal kaons, this selection has an efficiency of 
about 85% and a pion misidentification rate of less than 
2%, as determined from control samples of D* — > D°ir, 
D° — > Kit events. The detailed performance of the kaon 
selection has been characterized as a function of labora- 
tory momentum and can be seen in Fig. 4. 



We reconstruct the r\ in two final states: r\ — > 77 (?7 77 ) 
r~7r (77377)- For the rf , we reconstruct two 
7/ -> p° 7 (77^) and 77' -> i]tt + tt- (r^J, 
77 (except in the T]', 3 ^„„K*° mode, where we 



and 77 

final states: 
with 77 

also include 77 — > /t+/t~/t°). In the B° — > urn channel, we 
reconstruct ui — > 7t + /t~/t°; for _B° — > </>7r° we reconstruct 
<fi — > K + K~ . We place the following requirements on the 
invariant masses of the resonance candidates (in MeV): 
520 < m V37r < 570, 490 < m, q _ n < 600, 910 < Try < 1000 
for rj' pj and 77^, 735 < m w < 835, and 990 < < 
1050. These ranges can be seen graphically in Fig. 8 in 
Sec. VI B. The mass requirements for these resonances 
are loose to keep appropriate sidebands for fitting; the 
resonance shapes used for fitting are discussed in Sec. VI. 

For 77—8-77 candidates we require \Tlri\ to be less than 
0.86, where H v is the cosine of the 77 decay angle. The 
decay angle is defined, in the 77 rest frame, as the an- 
gle between one of the photons and the direction of the 
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FIG. 4: Identification (ID) efficiency of the charged kaon se- 
lection as a function of the kaon laboratory momentum P l £ h 
(top), and fraction of charged pions misidentified (mis-ID) as 
kaons as a function of the pion laboratory momentum P^ b 
(bottom). The error bars represent statisticai uncertainties 
in the control sample of kaons and pions from D* — > D°tt, 
D° — > Kit decays. 
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boost needed to get to this frame from the B center- 
of-mass (CM) frame. This requirement removes very 
asymmetric decays of the 77, where one photon carries 
most of the particle's energy. It is effective against high- 
energy background photons from B — ► K*-f that combine 
with a random low-energy photon to form an invariant 
mass in the range chosen for the 77 — > 77 decay. For the 
rj' n p + channel, the r\ — > 77 mass range is tightened to 
510 < m 77 < 580 MeV to reduce the continuum back- 
ground in the sample. 



C. Photon and 7r° selection 

Photons are reconstructed from energy depositions in 
the electromagnetic calorimeter which are not associated 
with a charged track. We require that all photon can- 
didates have an energy greater than 30 MeV except for 
the modes rf'^-K , ujtt , and (fin , where there is signifi- 
cant combinatorial background arising from low-energy 
photons. For these modes, we tighten the photon-energy 
requirement to 50 MeV for all photons. For r\ — > 77, we 
require each photon energy to be greater than 100 MeV, 
and for the rj — > p°j modes, we require the photon from 
the 77' decay to exceed 200 MeV. 

We select neutral-pion candidates from two photon 
clusters with the requirement that the 77 invariant mass 
satisfy 120 < m,« < 150 MeV. The mass of a ir° can- 
didate meeting this criterion is then constrained to the 
nominal value [22] and, when combined with other tracks 
or neutrals to form a B candidate, to originate from the 
B candidate vertex. This procedure improves the mass 
and energy resolution of the parent particle. 

For the primary tt° in r/'V decays, photon candidates 
are required to be consistent with the expected lateral 
shower shape, and the magnitude of the cosine of the n° 
decay angle (defined as for the 77) must be less than 0.95. 



D. K° selection 

For decay chains containing a K°, we reconstruct 
only the K® — > n + ir~ decay. The invariant mass of 
the candidate K® is required to lie within the range 
488 < m n + v - < 508 MeV. We also perform a vertex- 
constrained fit to require that the two tracks originate 
from a common vertex, and require that the lifetime 
significance of the {t/<j t ) be > 3, where a T is the 
uncertainty in the lifetime determined from the vertex- 
constrained fit. 



E. K* and p selection 



-*+ 



or 



We reconstruct the K* + as either K + tt° (K K+ 
K° s w+ (K*+ >7r+ ), and the K*° as K+tt- {K$ +lt _). The 
p + is reconstructed as tt+tt and the p° as tt+tt - . A 
vertex fit is performed when reconstructing the resonant 



K* or p candidate. We require the invariant masses (in 
MeV) of the resonance candidates to be in the ranges: 
755 < m K7T < 1035, 470 < m T+7r o < 1070, and 510 < 
m w + w - < 1060. The lower limit on the p a candidate 
invariant mass is chosen to reject background from K® 
decays. 

For decay chains involving a charged K* or p, we de- 
fine H, the cosine of the angle between the pion and the 
negative of the B momentum in the vector-meson rest 
frame. For p + decays, the direction is that of the tt°. 
For p° decays, we use only the magnitude of H, which 
is independent of the choice of reference pion. For these 
decays with a 7r° in the final state, we require that Ji be 
greater than —0.5 to reject combinatorial background. 



F. B meson selection 

A £>-mcson candidate is characterized kincmatically 
by the energy-substituted mass tties and by the energy 
difference AE, defined as 



m E s 



is + p • p B 



En 



and 



AE = (2q q B - s) /2y/7s , 



(2) 
(3) 



where qs = (Eb,Pb) and qo = (En, po) are the four vec- 
tors of the -B-candidate and the initial electron-positron 
system, respectively, and s is the square of the invariant 
mass of the electron-positron system. When expressed in 
the T(4S) frame, these quantities take the simpler but 
equivalent form 



AE 



1 *2 



and 



E* B 



(4) 
(5) 



where the asterisk denotes the value in the T(4S) frame. 
The mode-dependent resolutions on these quantities for 
signal events are about 3 MeV for tties, and 30-60 MeV 
for AE. 

We require 5.20 < m ES < 5.29 GeV and \AE\ < 0.2 
GeV for all but the r/^n , ujtt , and 4>-k° modes, where 
we loosen the AE range to \AE\ < 0.3 GeV to account 
for poorer detector resolution in these channels. 

When multiple B candidates from the same event pass 
the selection requirements, we choose a single candidate 
based on criteria described below. The average number 
of candidates per event depends on the mode; it is typ- 
ically about 1.2 and is always less than 1.5. We find 
that 70-90% of the events have a single combination and 
about 90% of the rest have two combinations. In decays 
containing an 77 and a K* or p, we select the candidate 
with the smallest \ 2 formed from the 77 and K* or p 
masses. For decays containing 77' — > rjir+TT^, the x 2 is 
formed from the masses of the 77' and 77 candidates. For 
all other decays, we retain the candidate that has the 
mass of the primary resonance (rj^\ lu, or <f>) closest to 
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the nominal value [22] . We have checked that this choice 
introduces no significant yield bias, in part because, for 
the primary resonance mass, there is an adjustable peak- 
ing component included in the fit, which would account 
for any small distortion due to this selection. 

IV. SOURCES OF BACKGROUND AND 
SUPPRESSION TECHNIQUES 

Production of BB pairs accounts for a relatively small 
fraction of the e + e~ cross section even at the peak of the 
T[AS) resonance. Upsilon production amounts to about 
25% of the total hadronic cross section, while tau-pair 
production and other QED processes occur as well. We 
describe below several sources of background, and discuss 
techniques for distinguishing them from signal. 

A. QED and tau-pair backgrounds 

Two-photon processes, Bhabha scattering, muon-pair 
production and tau pair production are characterized by 
low charged track multiplicities. Bhabha and muon-pair 
events are significantly prescaled at the trigger level. We 
further suppress these and other tau and QED processes 
via a minimum requirement on the event track multiplic- 
ity. We require the event to contain at least one track 
more than the topology of our final state, or three tracks, 
whichever is larger. We also place a requirement on the 
ratio of the second to the zeroth Fox- Wolfram moments 
[23], i?2 < 0.98, calculated with both charged tracks and 
neutral energy depositions. These selection criteria are 
more than 90% efficient when applied to signal. From 
MC simulations we have determined that the remaining 
background from these sources is negligible. 

B. QCD continuum backgrounds 

The primary source of background to all charmless 
hadronic decays of the B meson arises from continuum 
quark-antiquark production. The fact that these events 
are produced well above threshold provides the means by 
which they can be rejected, as the hadronization prod- 
ucts are produced in a jet-like topology. In strong con- 
trast, B mesons resulting from T(4S) decays are pro- 
duced just above threshold. Thus the final-state particles 
in the signal are distributed approximately isotropically 
in the CM frame. 

Several event-shape variables are designed to take ad- 
vantage of this difference. We define the thrust axis for 
a collection of particles as the axis that maximizes the 
sum of the magnitudes of the longitudinal momenta with 
respect to the axis. The angle 9t between the thrust 
axis of the B candidate and that of the rest of the tracks 
and neutral clusters in the event, calculated in the T(45) 
frame, is the most powerful of the shape variables we 
employ. The distribution of the magnitude of cos#t is 



sharply peaked near 1 for combinations drawn from jet- 
like qq pairs and is nearly uniform for the isotropic B- 
meson decays. This behavior is shown in Fig. 5. The 
selection criterion placed on cos 9t is optimized for each 
channel to maximize our sensitivity to signal in the pres- 
ence of continuum background and to reduce the size of 
the sample entering the fit. The optimization procedure 
is described in Sec. VII. The maximum allowed value of 
| cos#t| chosen for each signal mode is listed in Table I. 
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FIG. 5: Distribution in | cos 0t | for a typical B meson decay 
(B° — > T]' p7 7T° MC, solid points) and for the corresponding 
continuum background data (open circles). 

Further use of the event topology is made via the con- 
struction of a Fisher discriminant T ', which is subse- 
quently used as a discriminating variable in the likeli- 
hood fit. The Fisher discriminant we use is an optimized 
linear combination of the remaining event shape infor- 
mation (excluding cos 6>t , which we have already used in 
our preselection requirements). The variables entering 
the Fisher discriminant are the angles with respect to 
the beam axis of the B momentum and B thrust axis 
(in the T(4S) frame), and the zeroth and second angu- 
lar moments Lo,2 of the energy flow about the B thrust 
axis. The moments are defined by Lj = J2iPi x |cos0i| J , 
where 6i is the angle with respect to the B thrust axis of 
track or neutral cluster i, pi is its momentum, and the 
sum excludes the B candidate. The coefficients used to 
combine these variables are chosen to maximize the sep- 
aration (difference of means divided by quadrature sum 
of errors) between the signal and continuum background 
distributions of Lj, and are determined from studies of 
signal MC and off-peak data. We have studied the op- 
timization of T for a variety of signal modes, and find 
that the optimal sets of coefficients are nearly identical 
for all. Thus we do not re-optimize the Fisher coefficients 
for each individual decay. Because the information con- 
tained in T is correlated with | cos 6>t | , the separation be- 
tween signal and background is dependent on the | cos 9t \ 
requirement made prior to the formation of T . In Fig. 6, 
we show the Fisher-discriminant distribution for signal 
and continuum background for the B~ — > D°tt~ control 
sample. 
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FIG. 6: Distributions of Fisher-discriminant output for the 
data control mode B~ — > D°ty", D° — > K~ir + Tr a (points 
with error bars), corresponding signal Monte Carlo (solid his- 
togram), continuum data (open circles) and continuum Monte 
Carlo (dashed histogram) after requiring | cos 6t \ < 0.9. The 
Fisher discriminant and | cos 8t | are strongly correlated, so 
the separation depends on this requirement. 



C. BB backgrounds 



Most charmless hadronic-i?-decay analyses do not have 
much background from other B decays. Specifically, 
since most B mesons decay via b — > c transitions, the 
strange and light meson decay products from such de- 
cays result from b — > c — > q cascades, and thus have 
lower momentum than those expected in the signal final 
states. This small background is included in our qq back- 
ground PDF shapes (see next section) since the shapes 
are extracted from on-peak data. 

We have found, however, that some of the signal modes 
(see Table II in Sec. IX) do suffer from backgrounds from 
charmless hadronic decay modes. We investigate back- 
grounds that may not be completely suppressed by the 
selection criteria defined in Sec. Ill with Monte Carlo 
samples of BB events corresponding to several times the 
number of such events in the dataset. When we find 
an indication of a high selection rate for a particular 
background decay mode, we use the experimentally mea- 
sured (when available) or theoretically predicted branch- 
ing fraction of that mode to determine its expected con- 
tribution. Fits to simulated experiments such as those 
described in Sec. VII are used to evaluate whether such 
events cause a significant bias to the measured signal 
yield. Based on these studies, we have adjusted (while 
still blind) some selection criteria and in some cases 
added a component to the ML fit to account explicitly for 
the remaining BB background contributions. System- 
atic errors account for the uncertainties in this method. 
The details of this procedure are described below. 



V. MAXIMUM LIKELIHOOD FIT 

We use an unbinncd, extended maximum likelihood 
fit to extract signal yields for our modes. A subsample 
of events to fit for each decay channel is selected as de- 
scribed in Sec. III. The sample sizes for the decay chains 
reported here range from 700 to 30,000 events, where we 
include sidebands in all discriminating variables in order 
to parameterize the backgrounds. 

A. Likelihood function 

The likelihood function incorporates several discrimi- 
nating variables to distinguish signal from the large num- 
ber of background events retained by the sample selec- 
tion. We describe the -B-decay kinematics with two vari- 
ables: AE and ttzes (as defined in Sec. IIIF). We also 
include the mass of the primary resonance candidate (ro^, 
m v i, mK*, m p , m w , or m^) and the Fisher discriminant 
T . For the vector-pseudoscalar modes with a K* , p, to, 
or cf>, we also include in the fit the hclicity cosine H of 
the vector meson. For the K*, p, and </>, TL is defined in 
Sec. HIE. For the B° — > uitt decay, TL is defined as the 
cosine of the angle between the normal to the lo decay 
plane (the plane of the three pions in the lo rest frame) 
and the flight direction of the uj, measured in the lo rest 
frame. 

Because correlations among the discriminating vari- 
ables (except resonance mass and TL for background) 
in the selected data are small, we take the probabil- 
ity distribution function (PDF) for each event i to be 
a product of the PDFs for the separate discriminating 
variables. We define hypotheses j, where j can be sig- 
nal, continuum background, or (where appropriate) BB 
background. The PDFs can be written as 

V)=V 3 {^ ES )V J {^)V J {r)V 3 {r^ P )V {m\,,rC), (6) 

where mp indicates the pseudoscalar candidate mass in 
the fit (absent for B° — > uitt and B° — > 07r° modes) and 
my indicates the vector candidate mass (absent for the 
r/^ir modes). 

The likelihood function for each decay mode is 

where Yj is the yield of events for hypothesis j (to be 
found by the fitter) and N is the observed number of 
events in the sample. The first factor takes into account 
the Poisson fluctuations in the total number of events. 



VI. SIGNAL AND BACKGROUND MODEL 

We determine the PDFs for signal from MC distribu- 
tions for each discriminating variable. The PDFs for BB 
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background (where appropriate) arise from fitting the 
composite BB MC sample, described in Sec. VIA. For 
the continuum background we establish the functional 
forms and initial parameter values of the PDFs with data 
from sidebands in toes or AE. We then refine the main 
background parameters (excluding resonance-mass cen- 
tral values and widths) by allowing them to float in the 
final fit so that they are determined by the full data sam- 
ple. The following sections describe first the construction 
of samples to represent BB background, and then the 
control samples used to validate the PDF shapes and 
make adjustments to the means and widths of the distri- 
butions where needed. Finally we describe the detailed 
functional forms used to parameterize all of the signal 
and background distributions. 



A. Inclusion of BB background in the fits 

As discussed in Sec. IV C, backgrounds from other 
charmless B decays need to be accounted for explicitly 
in the maximum likelihood fit for some decay chains. 

Since we find that the signal yield bias due to BB 
background for the rj^K* channels is less than 1% of the 
signal yield, we do not include a BB component for these 
modes. For all modes with a K* + — > K + tt decay, nearly 
all BB backgrounds are removed by the requirement 7i > 
—0.5. This requirement is also helpful in reducing the 
BB background for decays with a p + — ► 7r + 7r°, though 
sufficient background remains to be included in the fit. 
For all other modes except B° — > (f>ir°, we include a BB 
component in the fit. The fit number of BB events is 
a small fraction of the total sample and is tabulated in 
Table II in Sec. IX. 

The PDFs for BB background are determined by fit- 
ting a sample of MC events composed of several charm- 
less decay chains, with the PDF shapes described be- 
low. For the r\ — > 77 channels, the BB background is 
dominated by B — ► K*j decays, even after the rj de- 
cay angle requirement, due to the relatively large if* 7 
branching fraction (40 x 10~ 6 ). For the B — ► r\p chan- 
nels, the largest backgrounds are from r]K* decays, with 
misidentification of the charged kaon or loss of the kaon 
while selecting a pion from the other B. For the rf chan- 
nels, the dominant backgrounds in all modes, except for 
Tj pl Pi arise from B — > rj'K decays, due to the relatively 
large branching fraction (~ 70 x 10~ 6 ). Another import- 
ant background for the rfp^K* channels, is K*p° decays, 
where the p is combined with a photon to fake an rf . For 
the r)' pl p and r]' p7 Tr modes, BB backgrounds are primar- 
ily from B + — > p + p° and B° — > p + p~ decays. For the 
decays with a primary ir°, the largest backgrounds are 
from B + — > rf 1 ' p + and B + — > up + decays, where due to 
the forward-backward peaking of the p + Ti. distribution, 
the 7T° is often energetic and the charged pion is lost. 



B. PDF corrections from data control samples 
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FIG. 7: Distributions of (a) toes and (b) AE from the B~ —> 
■k~ D° data sample used to determine the small corrections to 
signal Monte Carlo PDF shapes. 



We validate the simulation on which we rely for signal 
PDFs by comparing critical distributions of discriminat- 
ing variables in MC with those from large data control 
samples. For toes and AE (see Fig. 7), we use the decays 
B~ -> tt-D° and B~ -> p~ D° with D° -> A'-tt+tt , 
which have similar topology to the modes under study 
here. We select these samples by making loose require- 
ments on toes and AE, and more stringent selections on 
cos #t and the D° and p candidate masses (as appropri- 
ate). We also place kinematic requirements on the D and 
B daughters to force the charmed decay to look as much 
like that of a charmless decay as possible without elimi- 
nating the control-sample signal. These selection criteria 
are applied both to the data and to a MC mixture of re- 
lated B — ► DX and B — > D*X decays, which simulates 
the crossfeed from D* — » D° decays observed in data. 
From these control samples, we determine small adjust- 
ments to the mean value of the signal toes distribution 
and to the resolution of the AE distribution compared 
with Monte Carlo. For T we use parameters found from a 
sample of approximately 500 B + — > rj' pi K + events, with 
a cos 9t requirement matching that used for each signal 
mode. 

For the mass shapes of the resonances, we study inclu- 
sive resonance production in the off-peak data and corre- 
sponding continuum MC. In each sample, we reconstruct 
resonance candidates involved in our final states, requir- 
ing a minimum value of the candidate CM momentum 
of 1.9 GeV to reflect the kinematics of our final states. 
The resolutions and means of the invariant mass distribu- 
tions are compared, and we adjust the means and widths 
of PDF parameterizations based on the outcome of these 
results. A typical mass distribution for each resonance is 
shown in Fig. 8. 



13 



L (a) 












: 



0.5 0.52 0.54 0.56 0.58 0.6 
m„ [GeV] 




0.8 0.85 0.9 0.95 1 

m K . [GeV] 




1 (b) 

\ J 




0.92 0.94 


0.96 0.98 

m„. [GeV] 


; (d) 






0.6 


0.8 1 

m p [GeV] 


: (f) f 







1.03 1.04 
m t [GeV] 



FIG. 8: Distributions of the candidate masses for resonant 
decays from the on-peak sideband samples in data that are 
used to describe the signal PDF shapes (see Sees. VI F and 
VI G). For each distribution a real resonance signal compo- 
nent is evident above a combinatorial background component: 
(a and b) the four rf'> candidate mass combinations from the 
r/p and rj'p samples; (c) K* candidate mass from the T]'„„K* 
sample; (d) primary p candidate mass from the rj' pj p sample; 
(e) uj candidate mass from the lotv sample; (f) <f> candidate 
mass from the 4>tv° sample. In (a) the arrows indicate the 
narrower mass requirement for the r\ — > 7r + 7r~7r° decay. The 
same range is used even for the narrower vj — ► t]it + -k~ dis- 
tribution, shown as the lower plots in (b). For the K* and 
p cases, we do not show both charges since the distributions 
are very similar. 



C. ttxes parameterization 

The signal distribution in mES is parameterized by two 
Gaussian functions centered near the mass of the B me- 
son. The second Gaussian typically accounts for less than 
20% of the total area, and has a larger width to take into 
account the tails of the distribution, which arise primar- 
ily from misreconstructed signal events. In continuum 
background, we model wes by a phase-space-motivated 
empirical function [21] of the form 



f(x) oc x\ 1 — x 2 



exp 



-£(i-z 2 )], 



(8) 



where we define x = Im^sl yfs, an d £ is a parameter de- 
termined by the fit. In BB background samples, we find 
that the toes distribution is well-described by adding a 
simple Gaussian function to the empirical shape in Eq. 8; 



a similar alternate form of a Gaussian convolved with an 
exponential is used for some channels. 



D. A.E parameterization 

For AE, we fit the signal distribution with two Gaus- 
sian functions, both centered near zero. The broad Gaus- 
sian has a width about five times larger than the narrow 
Gaussian; this accounts for energy loss before or leakage 
out of the EMC, as well as incorrect candidate combina- 
tions in true signal events. The broad Gaussian compo- 
nent becomes larger as more of the final state energy is 
carried by neutral particles. The primary Gaussian func- 
tion accounts for about (60-80)% of the total area in all 
modes except rf''p where it is between 30% and 60%. 
For continuum background, we model the AE distribu- 
tion with a linear or quadratic polynomial as required by 
the data. The BB background is described well by two 
Gaussian functions peaking at negative (positive) AE, 
accounting for backgrounds that have a larger (smaller) 
number of tracks and neutrals in the final state than the 
signal. 



E. Fisher parameterization 

For both signal and background, the Fisher distribu- 
tion T is described well by a Gaussian function with 
different widths to the left and right of the mean. For 
the continuum background distribution, we also include 
a second Gaussian function with a larger width to ac- 
count for a small tail in the signal T region. This ad- 
ditional component of the PDF is important, because it 
prevents the background probability from becoming in- 
finitesimally small in the region where signal lies. As 
shown in Fig. 6, the mean of the continuum background 
distribution is approximately 2a greater than the mean 
of the signal peak, allowing for strong discrimination be- 
tween the two. Because T describes the overall shape 
of the event, the distribution for BB background is very 
similar to the signal distribution; hence this variable has 
little discriminating power against BB background. 



F. Pseudoscalar mass parameterization 

The pseudoscalar candidate mass distributions for sig- 
nal are described well by the sum of two Gaussian func- 
tions. We use MC values for the means and widths of 
these Gaussians, corrected where necessary by using sam- 
ples such as those shown in Fig. 8. In continuum back- 
ground, we fit the data with two Gaussian functions, 
where we fix the means and widths to those used for 
signal, and include a linear or quadratic term to account 
for non-resonant background. The fraction of resonant to 
non-resonant background is allowed to float in the final 
fit. When there is no discernible resonant component, 
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as in t]' — > p 7, floating this parameter can cause conver- 
gence issues in the final ML fit. If validation studies show 
this effect, the resonant fraction is fixed in the final anal- 
ysis. For BB background, wc use the same functional 
form as in continuum background; whether or not there 
is a true resonant component in BB background depends 
upon the charmless decay chains expected to contribute. 

G. Vector mass and helicity parameterization 

In pseudoscalar-vector decays of the B meson, the vec- 
tor meson has a helicity-angle distribution proportional 
to TL 2 for true signal events. We model the vector-meson 
helicity distribution for signal with a polynomial times 
a threshold function that allows for the effects of accep- 
tance. The signal K* and ui invariant-mass distributions 
are described by Breit-Wigner shapes. The <j> and p line 
shapes are found to be modeled well by two Gaussian 
functions; these do not fit well to a Breit-Wigner shape 
because of non-negligible mass resolution ((f)) or misre- 
constructed p candidates in real signal events (p) . For the 
p and other wide distributions there is as much as 10% 
loss of efficiency due to the effect of the mass range re- 
quirements; this effect is included in the overall efficiency 
estimate and its uncertainty is included in systematic er- 
rors discussed in Sec. X. See Fig. 8 for illustrations of 
these distributions. 

Because the shape of the helicity angle can be differ- 
ent for continuum background with and without a true 
vector resonance, we use a two-dimensional PDF to de- 
scribe the resonance mass distribution and the helicity- 
angle distribution. We would expect that the background 
H would have a nearly uniform distribution, correspond- 
ing to a sum of combinatorial resonance background and 
background of true resonances from various production 
mechanisms. We find that the pure-background shape is 
modeled well by a second order polynomial with only a 
small amount of curvature and the true-resonance com- 
ponent is a separate low-order-polynomial shape. The 
mass parameters for the true-resonance component are 
fixed to be the same as for the signal. 

The BB background component of H is modeled by 
a single fourth-degree polynomial. We parameterize the 
resonance mass distribution with two Gaussian functions 
plus a linear or quadratic polynomial, allowing the means 
and widths of the Gaussians to float if the resonant com- 
ponent of the background differs from the signal reso- 
nance. This is especially necessary when BB background 
arises when a misidentified kaon from a K* causes its re- 
construction as a p. 

The requirement that charged tracks have pr > 100 
MeV (Sec. Ill) can induce a "roll-off" effect near TC val- 
ues of ±1. In particular, for decays of a K* or p with 
a charged pion, the helicity distribution of the vector 
meson shows a characteristic roll-off in the region popu- 
lated by low-momentum pions. This effect is absent for 
charged kaons since there are no kaons with px < 100 
MeV. We model the roll-off in both the signal and back- 



ground Ti. distributions by multiplying the primary PDF 
shape by an appropriate Fermi-Dirac threshold function. 
The parameters of this roll-off function are constrained to 
be the same for signal and both background components. 
Because the ui helicity angle is defined from a three-body 
decay (u) — ► 7r + 7r~7r ), there is little correlation between 
low-momentum pions and helicity angle, and hence no 
significant roll-off. 



VII. FIT VALIDATION 

Before applying the fitting procedure to the data to 
extract the signal yields we subject it to several tests. 
Internal consistency is checked with fits to ensembles of 
"experiments" generated by Monte Carlo from the PDFs. 
From these we establish the number of parameters associ- 
ated with the PDF shapes that can be left free in addition 
to the yields. Ensemble distributions of the fitted pa- 
rameters verify that the generated values are reproduced 
with the expected resolution. The ensemble distribution 
of In C itself provides a reference to check the goodness of 
fit of the final measurement once it has been performed. 

We account for possible biases due to neglecting cor- 
relations among discriminating variables in the PDFs by 
fitting ensembles of experiments into which we have em- 
bedded the expected number of signal events randomly 
extracted from the detailed MC samples, where correla- 
tions are modeled fully. We find a positive bias of a few 
events for most modes, as shown in Table I. Events from 
a weighted mixture of simulated BB background decays 
are included where significant, and so the bias we mea- 
sure includes the effect of crossfeed from these modes. 

For modes with low background and small signal 
yields, the ensemble yield distribution may exhibit a sig- 
nificant negative tail. This is due to the nature of the 
maximum likelihood method, which is known to be bi- 
ased for small samples. The source of the bias is the in- 
sufficient number of events for which the probability for 
the signal hypothesis is larger than the probability for the 
background hypothesis. This results in a negative bias, 
which is taken as the mean of the yield distribution from 
the fits to the ensembles described above. Examples of 
modes with negative bias can be found in Table I. By 
subtracting the bias we correct for this effect on average, 
and we include the uncertainty as a systematic error. 

This same procedure for generating and fitting sim- 
plified MC samples is used to find an optimal selection 
requirement for the cos 6> T variable in the early stages of 
each analysis. The studies are performed for a range of 
selection values, to minimize the fractional error on the 
signal yield. The optimal values of the cos6>t require- 
ment that are chosen are given in Table I. 

Finally, we apply the fit to the off-peak data to confirm 
that we find no fake signals in a sample with no signal 
events. 
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TABLE I: For each B decay chain we present the optimized 
| cos # T | requirement, the number of on-peak events passing 
the preselection requirements, and the fit bias Yj, determined 
from simulated experiments (the uncertainty on this bias is 
discussed in Sec. X). 



Mode 


Max | cos 9t \ 


#Events Fit bias, Y b 






in fit 


(events) 


B+ -> r,K*+ 










0.90 


7573 


4.7 




0.90 


4132 


1.7 




0.90 


4Q74 


1 




0.90 


2835 


0.3 


B° -» n K* 1 










0.90 


12179 


8.1 




0.90 


6440 


1.8 


B+ 7JP+ 










0.80 


17084 


1.3 




0.90 


16106 


1.0 


B° -> V p° 










0.70 


11107 


-1.0 


m^p° 


0.80 


8347 


2.3 


B° -» V n° 








»7 77 7r 


0.80 


5379 


-1.1 


7737T7T 


0.80 


2271 


0.7 



B+ -> riK*+ 








ri K* + 
ri K*^ 


0.90 


2973 


-4.5 


0.75 


13299 


3.6 


ri K* + 


0.90 


2009 


0.0 


ri K* + 


0.75 


8205 


0.6 


b° -» riK*° 








ri K*° 


0.90 


4808 


-3.7 


t IS*® 


0.90 


695 


1.7 


0.75 


20504 


4.2 


B+ -» rip+ 








1 + 


0.90 


8737 


2.1 


Vp~ f P + 
B° -+ rip" 


0.65 


28933 


7.8 


0.90 


9515 


-3.7 


B° -» V 'tv° 










0.90 


3491 


-3.5 


»7p 7 7T° 


0.70 


11426 


2.8 



B° - oj-k 


0.80 


18986 


-2.1 


B° -» 4m° 


0.90 


4840 


-1.1 



VIII. EFFICIENCIES AND EFFICIENCY 
CORRECTIONS 

The efficiency is determined by the ratio of the number 
of signal Monte Carlo events passing preselection to the 
total number of generated MC signal events. This effi- 
ciency is corrected for differences between the true detec- 
tor efficiencies and those simulated in Monte Carlo. From 
a study of absolute tracking efficiency, we apply a correc- 
tion of (1-7)%, depending on the number of charged par- 
ticles in the decay channel and assign a systematic error 
of 0.8% per track. The Kg efficiency correction is taken 



from an independent study of the vertex-displacement 
dependence of the efficiency for inclusive samples of K® 
mesons from the data and from MC. The overall cor- 
rection for the topologies represented by our decays is 
0.971 ± 0.030. For the six decays with a primary tt° and 
the four with a K* + or p + decaying to a final state with 
an energetic ir°, we determine a correction from a sample 
of tau decays. For these cases, the n° efficiency is 6-11% 
lower for data than MC. 



IX. FIT RESULTS 

The branching fraction for each decay chain is obtained 
from 



B 



Y-Y b 



(9) 



where Y is the yield of signal events from the fit, Y, is 
the fit bias discussed in Sec. VII and given in Table I, 
e is the efficiency, Bi is the branching fraction for the 
ith unstable B daughter (Bi having been set to unity in 
the MC simulation), and Nb is the number of produced 
B + or B° mesons. The values of Bi are taken from Par- 
ticle Data Group world averages [22]. The number of 
produced B mesons is computed with the assumption of 
equal production rates of charged and neutral B pairs 
[24]. 

In Table II, we show the results of the final ML fits to 
the on-peak data, with the yields for signal and BB back- 
ground, where applicable. The latter is often uncertain 
due to the large correlation with the qq background com- 
ponent, but this uncertainty is not problematic because 
the correlation with signal is small. We also show the effi- 
ciencies, daughter branching-fraction products, and esti- 
mated effective purity of the sample. We report the sta- 
tistical significance for the individual decay chains and 
display the significance including systematic uncertain- 
ties for the combined result in each channel. The purity 
is the ratio of the signal yield to the effective background 
plus signal; we estimate the denominator by taking the 
square of the uncertainty of the signal yield as the sum of 
effective background plus signal. Where the signal yields 
are small the purity is not very meaningful, so we do not 
report the purity if it is below 10%. Branching fractions 
are given for individual fits to each submode as well as 
the result of combining several submodcs. Since the lat- 
ter procedure involves systematic as well as statistical 
errors, we defer the description to Sec. XL The final col- 
umn in Table II gives the charge asymmetry (A c h), as 
defined in Sec. I. 

The statistical error on the yield is given by the change 
in the central value when the quantity — 21n£ increases 
by one unit. The statistical significance is taken as the 
square root of the difference between the value of —2 In £ 
for zero signal and the value at its minimum. The 90% 
C.L. upper limit quoted in Sec. XIII is the solution 690 
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TABLE II: Fitted event yields (Y = signal yield, Y BB - = BB yield), purity (P, see text), efficiency (e), daughter product 
branching fractions (in percent), significance S(a) (which includes systematic errors), fit branching fractions, 90% C.L. upper 
limits, and charge asymmetries. Also shown are the results of combining daughter decay chains where more than one contribute. 
For the final branching fraction and charge asymmetry results, the systematic errors are also given. 



Mode Y Y bE P(%) e T\Bi S{a) B(1(T 6 ) UL (f(T 6 ) A ch 



B+ -f r]K*+ 9 25.6 ±4.0 ±2.4 +0.13 ± 0.14 ± 0.02 





46 ± 12 


25 ± 15 


35 


24.0 


9.0 


4.9 


22 ± 6 




±0.03 ± 0.24 




27 ± 8 




45 


17.1 


5.2 


5.0 


33 ± 10 




+0.46I 028 




30 ± 9 




45 


8.8 


13.1 


5.7 


29 ± 8 




—0.11 ± 0.28 




10 ± 5 




43 


6.6 


7.5 


3.2 


22 ± 11 




+0 37+21? 


B° -> r)K* 












11 


18.6 ± 2.3 ± 1.2 




±0.02 ± 0.11 ± 0.02 




125 ± 16 


5 ±19 


50 


24.4 


26.3 


10.1 


20 ±3 




+0.12 ±0.13 




32 ±9 




47 


16.5 


15.1 


5.0 


14 ±4 




-0.39 ±0.25 


B+ -f ?7p + 












3.5 


9.2 ± 3.4 ± 1.0 


<14 




r7 77 p+ 


32 ±15 


-3 ±19 


14 


10.7 


39.4 


2.5 


8±4 








21 ± 11 


3 ± 11 


17 


8.6 


22.6 


2.4 


12 ±6 






B° - V p 












— 


-i.it2:S±o.4 


<1.5 




V~nP° 


-18 ± 18 


67 ± 38 


<10 


27.1 


39.4 




— 2 ± 2 






n 

m^p 


— 2 ± 4 


26 ± 10 


<10 


18.2 


22.6 




— 1 ± 1 






B° -> T)TT° 












0.8 


0.7+1:1 ± 0-3 


<2.5 






1±7 


-2 ±9 


<10 


19.3 


39.4 


0.3 


0± 1 








8 ± 7 


-8 ±5 


15 


14.9 


22.6 


1.1 


2±2 






B+ -f r)'K* + 












1.9 


6.3+t« ± 1.8 

— 3.6 


<14 




n' K* + 

17)7777 K®77 J <~ 

'P7 K n 7r + 

n' K""+ 


-8 ±4 


29 ±11 


<10 


17.5 


4.0 




-5 ±6 






16 ±9 


17 ± 12 


22 


13.5 


6.8 


1.7 


15 ± 11 






3±3 




13 


7.0 


5.8 


1.7 


8±7 






5±7 




<10 


5.6 


9.8 


0.6 


8 ±14 






B° _> rj'K* 












2.1 


4 1+ 21 ±12 


<7.6 




n' K*° 

'*7(77) 7r7r 

' 7^*0 


0±4 


18 ± 10 


<10 


17.8 


11.6 


1.0 


2±2 






I 1 _1_ c 

II ± O 


1 O _1_ A 

to ± y 


A *7 

4/ 


loo 
12. z 


6. ( 


o n 
2.U 


13 ± I 








15 ±10 


80 ±25 


17 


14.0 


19.7 


1.3 


5±4 






B+ ->■ T/V+ 












2.6 


12. 9+5 5 ± 2 -° 


<22 




'HrjiYiYp 


16 ±8 




25 


8.4 


17.5 


2.1 


11 ±6 






v'p~ f p + 


48 ±23 


61 ± 100 


<10 


6.5 


29.5 


1.7 


24 ± 13 






B° - r,'p° 


-1 ±4 


53 ±21 


<10 


19.7 


17.5 


0.5 


O.8+12 ± °- 9 


<4.3 
















0.7 


i.o±l:t ± °- 8 


<3.7 




7)' 7T° 

ITjlYTY " 


-2 ±3 


-8 ±4 


<10 


18.5 


17.5 


0.4 


1 ± 1 








17 ± 14 


-38 ± 78 


<10 


13.9 


29.5 


1.1 


4±4 






B° -> ^7T° 


-9 ±8 


9 ±18 


<10 


15.9 


89.1 




-0.6+Q 5 ± °- 2 


<1.2 




B° -> 07T° 


2±4 




<10 


28.6 


49.2 


0.7 


0.2^ ±0.1 


<1.0 





to the equation 



f °°£(b)db 



0.9, 



(10) 



where C(b) is the value of the maximum likelihood for 
branching fraction b. 

In Figs. 9-11 we show projections of all fit discriminat- 
ing variables for the rjK* and r)p + modes. Points with 
errors represent data, solid curves the full fit functions, 
and dashed curves the background functions. Since the 
r\ — > 77 and r\ — > 7r+7r~7r° components have very different 



resolutions, for the //-candidate mass plots we indicate 
with a dashed curve the full fit without the r\ — > 7r+7r~7r° 
signal component. We make these plots by selecting 
events with the ratio of signal to total likelihood (com- 
puted without the variable shown in the figure) exceeding 
a mode-dependent threshold that optimizes the expected 
sensitivity. The selection retains a fraction of the signal 
yield averaging about 70% across the decay sequences. 
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X. SYSTEMATIC UNCERTAINTIES 

We itemize estimates of the various sources of sys- 
tematic errors important for these measurements. Ta- 
bles III, IV, and V show the results of our evaluation of 
these uncertainties. We tabulate separately the additive 
and multiplicative uncertainties. That is, we distinguish 
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FIG. 9: Projections of the B-candidate discriminating vari- 
ables for B + — > rjK* + : (a) jties; (b) rj candidate mass; (c) 
AE; (d) K* candidate mass; (e) Fisher discriminant output; 
and (f) K* + helicity. See text for explanation of the points 
and curves. 




: (b) , 

t r 




1 ■ 1 ■ l 





5.2 5.22 5.24 5.26 5.28 
m ES (GeV) 

80F 



0.52 0.53 0.54 0.55 0.56 0.57 




0.1 0.2 
AE (GeV) 



0.9 1 
m K .(GeV) 




(0 
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-0.5 



0.5 1 
K* Helicity 



FIG. 10: Projections of the B candidate discriminating vari- 
ables for B° — > r/K* : (a) tties; (b) t] candidate mass; (c) 
AE; (d) K*° candidate mass; (e) Fisher discriminant output; 
and (f) K*° helicity. See text for explanation of the points 
and curves. 
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l*u , T T , — . — , 
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i.2 5.22 5.24 5.26 5.28 
m ES (GeV) 



0.52 0.53 0.54 0.55 0.56 0.57 
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FIG. 11: Projections of the B candidate discriminating vari- 
ables for B + — > rjp + : (a) thes; (b) r\ candidate mass; (c) AE; 
(d) p + candidate mass; (e) Fisher discriminant output; and 
(f) p + helicity. See text for explanation of the points and 
curves. 



those errors affecting the efficiency and total number of 
BB events from those that concern the bias of the yield, 
since only the latter affect the significance of the result. 
The two types of errors are comparable for modes with 
substantial yields but the additive errors dominate when 
the yields are small. Additionally we distinguish between 
those uncertainties that are correlated among different 
daughter decays of the same mode (C), and those that 
are uncorrelated (U). This distinction is relevant when 
multiple decay chains are combined (see Sec. XI). The 
final row of the table provides the total systematic error 
in the branching fraction for each of the submodes. 



A. Additive systematic errors 

Fit yield (U): Uncertainties due to imprecise knowledge 
of the background PDF parameters are included in the 
statistical errors since the main parameters are allowed 
to vary in the nominal fits. We have investigated the 
small correlations among background parameters and 
find these to have a negligible effect on signal yields. We 
include the uncertainty for the signal PDF parameters 
by determining the yield variations as individual param- 
eters are varied by uncertainties determined from fits to 
independent control samples (see Sec. VI B). 
Fit bias (U): This uncertainty is taken from the valida- 
tion procedure described in Sec. VII. We combine in 
quadrature terms, in order of relative importance, from 
(a) the positive bias (due to parameter correlations), (b) 
the negative bias for small event yields, (c) a small con- 
tribution from the modeling of the combinatorial compo- 
nent in signal, and (d) the statistical uncertainty in the 
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TABLE III: Estimates of systematic uncertainties in the branching fraction for B — > r\K* and B — + r\p decays. We distinguish 
between additive and multiplicative errors as well as errors that are correlated (C) or uncorrelated (U) among the submodes. 



(Quantity 




r]K*+ 








+ 




t] decay 


77 


Sty 


77 


Ovl 


77 


3tt 


77 


OTT 


77 


3tt 


K , p decay 


l\ TV 




K + 7Y° 


l\ 7T 


K + TY- 


K+-R- 


7T 71 


"7T + -77- 

7T 71 


7T + 7r" 




Additive errors (events) 






















Fit yield (U) 


2.1 


0.7 


0.4 


0.3 


2.2 


2.3 


2.3 


0.2 


5.2 


0.9 


Fit bias (U) 


2.5 


1.2 


0.3 


0.3 


4.2 


t.t 


1.1 


0.9 


0.3 


0. 1 


BB background (U) 


1.1 


0.3 


0.6 


0.5 


0.9 


0.3 


0.7 


0.4 


6.5 


1.2 


total additive (events) 


3.4 


1.4 


0.6 


0.5 


4.8 


2.6 


2.6 


1.0 


8.1 


1.6 


Multiplicative errors (%) 






















Tracking eff/qual (C) 


0.8 


2.4 


0.8 


2.4 


1.6 


3.2 


0.8 


2.4 


1.6 


3.2 


K° s efficiency (C) 


4.0 


4.0 


















TVfl Pk" mil l^iTiliPii"^/" I I; 1 

J_l cLL-lv lliLllLrlJJlH^lL^y 1 V_y J 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


7T / 7 eff (C) 


5.1 


5.1 


10.3 


10.3 


5.1 


5.1 


10.3 


10.3 


5.1 


5.1 


Number BB (C) 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Branching fractions (U) 


0.7 


1.8 


0.7 


1.8 


0.7 


1.8 


0.7 


1.8 


0.7 


1.8 


MC statistics (U) 


0.9 


1.1 


1.5 


1.8 


0.9 


1.1 


1.3 


1.7 


0.8 


1.0 


cos Ot (C) 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 


1.0 


0.5 


2.3 


1.0 


PID (C) 






1.0 


1.0 


1.0 


1.0 










Total multiplicative (%) 


6.8 


7.4 


10.6 


11.0 


5.8 


6.6 


10.6 


11.0 


6.1 


6.6 


Total a [B(10" 6 )] 


2.2 


2.9 


3.6 


2.6 


1.4 


1.5 


1.0 


1.5 


0.9 


0.4 



TABLE IV: Estimates of systematic uncertainties in the branching fraction for B — * rj'K* and B —> rj'p decays. The notation 
is the same as for Table III. 



Quantity 
r\ decay 
K* , p decay 


777T7T 
K°7Y + 


v' 

PI 


K*+ 

rjTTTV 
f K + 7Y 


PI 


■q-KTY 


7]'K*° 


PI 
K+n- 


v'p + 

TJ7V7T P7 
7T + 7T 7T+7T 


v'p° 

Tj-KTY 


Additive errors (events) 






















Fit yield (U) 


0.9 


1.3 


0.2 


0.8 


0.5 


2.8 


0.8 


0.8 


3.1 


1.0 


Fit bias (U) 


2.3 


1.9 


0.3 


0.4 


1.7 


0.8 


2.1 


1.2 


4.3 


1.9 


BB background (U) 


1.3 


0.9 


0.5 


0.5 


2.5 


0.5 


1.4 


1.2 


10.0 


1.8 


Total additive (events) 


2.8 


2.5 


0.6 


1.0 


3.1 


3.0 


2.6 


1.9 


11.3 


2.8 


Multiplicative errors (%) 






















Tracking eff/qual (C) 


2.4 


2.4 


2.4 


2.4 


3.2 


4.8 


3.2 


2.4 


2.4 


3.2 


K° s efficiency (C) 


4.0 


4.0 


















Track multiplicity (C) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


7T°/ 7 Off (C) 


5.4 


2.5 


10.4 


7.6 


5.4 


5.4 


2.5 


10.4 


7.6 


5.4 


Number BB (C) 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Branching fractions (U) 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


3.4 


MC statistics (U) 


1.0 


1.2 


1.7 


1.9 


1.1 


1.1 


1.1 


1.5 


1.8 


0.9 


COS0T (C) 


0.5 


1.8 


0.5 


1.8 


0.5 


1.4 


1.4 


0.5 


3.0 


0.5 


PID (C) 






1.0 


1.0 


1.0 


1.0 


1.0 








Total multiplicative (%) 


8.1 


6.8 


11.5 


9.2 


7.5 


8.4 


5.4 


11.4 


9.5 


7.4 


Total a [B(10" 6 )] 


4.5 


3.2 


1.9 


2.1 


1.7 


4.2 


1.1 


1.9 


6.9 


0.9 
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TABLE V: Estimates of systematic uncertainties in the 
branching fraction for the decays B° — ► r)ir , B° — > rj'ir , 
B° — > loty and B° — > 07r°. The notation is the same as for 
Table III. 



Quantity 




t/tt 


0JTV° 


07r° 


7/'' decay 


77 


3tt 


rj-Kir 


P7 






Add. err. (evts.) 














Fit yield (U) 


1.0 


1.3 


1.2 


2.8 


1.4 


0.8 


Fit bias (U) 


1.0 


0.3 


2.0 


1.3 


1.5 


0.7 


BB Bkg (U) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Total add. (evts.) 


1.7 


1.7 


2.5 


3.2 


2.3 


1.5 


Mult. err. (%) 














Track off. (C) 


0.0 


1.6 


1.6 


1.6 


1.6 


1.6 


Multiplicity (C) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


tt / 7 eff (C) 


14.9 


11.5 


13.1 


9.1 


11.7 


7.9 


Number BB (C) 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Br. frac. (U) 


0.8 


1.8 


3.4 


3.4 


0.8 


1.4 


MC stats. (U) 


1.2 


1.3 


1.2 


0.5 


1.2 


0.9 


COS0T (C) 


0.5 


0.5 


0.5 


1.3 


0.6 


0.5 


Total mult. (%) 


15.1 


11.9 


13.8 


10.1 


12.0 


8.4 


Total a [B(HT 6 )] 


0.3 


0.6 


0.9 


1.0 


0.2 


0.1 



determination of the bias. The first uncertainty (a) is 
taken to be one half of the positive bias, and the second 

(b) to be one half of the difference between the peak and 
mean yields of the ensemble distributions. Contribution 

(c) is small for all modes; we determine it using a com- 
parison of Monte Carlo and data for the B~ — > n~D° 
control sample. 

BB background (U): The BB background component, 
included in the fit for most decay chains, accounts for 
most uncertainties from BB background. We assign an 
additional uncertainty to account for modeling of this 
background. For the high-background rj' p + decay this 
involves explicit variation of the model. For the other 
modes it is taken to be 50% of the difference in the sig- 
nal yields when background is varied by its uncertainty 
(100% of the estimated effect when a BB background 
component is not included in the fit) and a contribution 
to account for uncertainty in the effect of the b — ► c back- 
ground. 



B. Multiplicative systematic errors 

Track finding/efficiency (C): As described in Sec. VIII, 
we assign a systematic error of 0.8% for each track (ex- 
cept for those from K® decays - see below). 
K® reconstruction efficiency (C): The K° efficiency sys- 
tematic uncertainty is taken from the study described in 
Sec. VIII with the addition of a contribution for recon- 
struction of the daughter charged tracks, giving a total 
uncertainty of 4% for decays with a K° in the final state. 
Track multiplicity (C): The inefficiency of the preselec- 



tion requirements for the number of tracks in the event 
is a few percent. We estimate an uncertainty of 1% from 
the uncertainty in the low-multiplicity tail of the B decay 
model. 

7, 7T°, r\~ n reconstruction efficiency (C): This uncertainty 
is estimated to be 2.5%/photon from a study of tau de- 
cays to modes with 7r°'s. For 7r°'s with energy greater 
than 1 GeV, there is an additional contribution to the 
uncertainty due to the overlap of the two showers, also 
evaluated from tau decays. 

Luminosity, B counting (C): From a sample of e + e~ — > 
fi + n~ decays, we estimate the uncertainty on the number 
of produced BB pairs to be 1.1%. 

Branching fractions of decay chain daughters (U): This 
is simply taken as the uncertainty on the daughter par- 
ticle branching fractions from Ref. [22]. 
MC statistics (U): The uncertainty due to finite signal 
MC sample sizes (typically 40,000 generated events) is 
given in the table. 

Event shape requirements (C): The uncertainties in the 
Fisher distribution T are included in the fit yield sys- 
tematic variation (see below). Uncertainties due to the 
cos6*t requirement are estimated to be one-half of the 
difference between the observed signal MC efficiency for 
the cos #t requirement used for each analysis and the 
expectation for a flat distribution. 

PIP (C): The uncertainties due to PID vetoes are negligi- 
ble. For analyses with a charged kaon, we estimate from 
independent samples an average efficiency uncertainty of 
1.0%. 



C. Charge asymmetry systematic errors 

For the B — ► r\K* analyses, the charged K used to 
define the asymmetry has a broad momentum spectrum. 
Auxiliary tracking studies place a stringent bound on de- 
tector charge-asymmetry effects at all momenta. Such 
tracking and PID systematic effects were studied in de- 
tail for the analysis of B — > (f)K* [25]. We assign the 
same 2% systematic uncertainty for A c h that was deter- 
mined in that study. In addition, we observe that the 
charge asymmetry of the continuum background is con- 
sistent with zero in all cases with a combined uncertainty 
below 1%. Finally we have measured the charge asym- 
metry for a control sample of B~ — ► D°p~ decays and 
find the result to be consistent with zero asymmetry, as 
expected. 



XI. COMBINED RESULTS 

To obtain the final results, we combine the branching 
fraction and charge asymmetry measurements from the 
individual daughter decay chains. The joint likelihood is 
given by the product, or equivalently — 2 ln£ is given by 
the sum, of contributions from the submodes. The sta- 
tistical contribution comes directly from the likelihood 
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fit, which reflects the non-Gaussian uncertainty associ- 
ated with small numbers of events. Before combining, 
we convolve each statistical C with a Gaussian function 
representing the part of the systematic error that is un- 
correlated among the submodes. The — 21n£ distribu- 
tions without systematic uncertainties give the combined 
statistical errors, while the distributions including corre- 
lated systematic uncertainties, give the total statistical 
and systematic errors. 

The resulting branching fractions and charge asym- 
metries are included in Table II, where the significance 
includes systematic uncertainties. 



XII. DISCUSSION 



More than six years have passed since the first report 
of a very large branching fraction for the decay B — ► rj'K, 
published in Ref. [2]. While it was expected [3] that the 
branching fraction for this decay and B — > r/K* would 
be relatively large and B — > r\K and B — ► rj'K* would 
be much smaller, most theoretical calculations could not 
account for a branching fraction as large as was mea- 
sured. The experimental situation with B — > vj ' K has 
remained largely the same even with quite precise new 
measurements; see for example Ref. [16]. The results 
presented in this paper complete the measurement of the 
four (r],r]')(K, K*) final states with a sensitivity in the 
branching fraction of a few times 10 -6 . The B — ► rjK* 
decays are found to have rather large branching fractions 
as expected and as first seen by CLEO [12]. BABAR has 
recently observed B + — ► r)K + for the first time [17] and 
finds the expected small branching fraction. We find no 
significant signal for B — > r]'K*, and the 90% C.L. up- 
per limit is not yet precise enough to determine whether a 
flavor-singlet component is present for this decay, though 
we do restrict the size of such a contribution. Such a sin- 
glet component (see Fig. Id) has been proposed as a par- 
tial explanation for the large rate for B — > rj'K by many 
authors, though with the restrictive limits for B — > rfK*, 
this now seems unlikely to play a significant role [26] . 

We also have evidence for the decay B + —> i]p + with a 
significance of 3.5a. We find no other significant signals 
and calculate upper limits for B + — > i]'p + and all of the 
neutral B decays with a p or ir° meson. This pattern is as 
expected since the penguin contribution in these decays 
is CKM suppressed and there is no external tree diagram 
for the B° decays. 

For the decays where we find significant signals, we 
also measure the charge asymmetry, which we find to 
be consistent with zero for all of these decays. These 
measurements are in agreement with the theoretical ex- 
pectations discussed in Sec. I and rule out substantial 
portions of the physical region. 



XIII. SUMMARY OF RESULTS 

We report measurements of branching fractions and 
charge asymmetries for _B-meson decays to r\ or rf with 
a K* , p, or 7T° as well as those channels with an u or <f> 
and a tt°. We find signals with high statistical signifi- 
cance in the B — > rjK* channels. We have evidence for 
the decay B + — > Tjp + (with significance 3.5a), which has 
not been seen previously. For branching fractions with 
significance less than four standard deviations, we quote 
both central values with errors and 90% C.L. upper lim- 
its. The observed values in the r\ channels are 

B{B + ^riK* + ) = (25.6 ± 4.0 ± 2.4 ) x 10" 



B(B° 



V K*°) = (18.6 ± 2.3 ± 1.2) x 10" 6 . 

B(B + ^rip + ) = (9.2 ± 3.4 ± 1.0) x 10~ 6 

< 14 x 10~ 6 , 
B(B Q ^rip°) = (-l.ll° j±0.4) x 10" 6 

< 1.5 x 10~ 6 , 

B(B° -> r/7r°) = {0.7±li ± 0.3) x 10~ 6 

< 2.5 x 10~ 6 . 



For the 77' channels, we find 

B(B + ^r)'K*+) = (6.3t£-.6 ± 1.8) x 10~ 6 

< 14 x 10~ 6 , 
B(B a ^ V 'K* Q ) = (4.1±?i ± 1.2) x 10" 6 

< 7.6 x 10" 6 , 

B(B+ -» 77> + ) = (12.9±|;§ ± 2.0) x 10~ 6 

< 22 x 10~ 6 , 

B(B° -» v 'p°) = {0.8t\i ± 0.9) x 10~ 6 

< 4.3 x 10~ 6 , 
B(B Q ^rj'7r°) = (l.Ot^ ± 0.8) x 10" 6 

< 3.7 x 10~ 6 . 

In the modes with a vector meson and a ir°, we observe 

B(B°^lutt°) = (-0.6I°;£ ± 0.2) x lO" 6 , 

< 1.2 x 10" 6 , 

B{B° -> </>7r°) = (0.2±g;| ±0.1) x lO" 6 , 

< 1.0 x 10" 6 . 

The results for B° — » uir° supersede the previous BABAR 
measurement of for this channel [15]. All of these results 
are substantially more precise than previous measure- 
ments from CLEO [12]. 

For the modes with significant signals, we measure the 
charge asymmetries 

A ch (r]K* + ) = +0.13 ±0.14 ±0.02, 
A ch (r]K* a ) = +0.02 ±0.11 ±0.02. 
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